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Abstract 

Membranous nephropathy (IVIN), characterized by the presence of diffuse thickening of the glomerular basement membrane 
and subepithelial in situ immune complex disposition, is the most common cause of Idiopathic nephrotic syndrome in adults, 
with an incidence of 5-10 per million per year. A number of studies have confirmed the relevance of several experimental 
insights to the pathogenesis of human MN, but the specific biomarkers of MN have not been fully elucidated. As a result, our 
knowledge of the alterations in histone methylation in MN is unclear. We used chromatin Immunoprecipitation followed by high- 
throughput sequencing (ChlP-seq) to analyze the variations in a methylated histone (H3K9me3) In peripheral blood 
mononuclear cells from 10 MN patients and 10 healthy subjects. There were 108 genes with significantly different expression In 
the MN patients compared with the normal controls. In MN patients, significantly increased activity was seen In 75 H3K9me3 
genes, and decreased activity was seen In 33, compared with healthy subjects. Five positive genes, DiGeorge syndrome 
critical region gene 6 (DGCR6), sorting nexin 16 (SNX16), contactin 4 (CNTN4), baculoviral lAP repeat containing 3 (BIRC3), 
and baculoviral lAP repeat containing 2 (BIRC2), were selected and quantified. There were alterations of H3K9me3 in MN 
patients. These may be candidates to help explain pathogenesis in MN patients. Such novel findings show that H3K9me3 may 
be a potential biomarker or promising target for epigenetic-based MN therapies. 
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Introduction 

Membranous nephropathy (MN), characterized by the 
presence of diffuse thickening of the glomerular basement 
membrane and subepithelial in situ immune complex 
disposition, shows a set of distinct morphologic features in 
glomeruli, such as glomerular basement membrane 
spikes, that form along the periphery of all glomerular 
capillary loops, and subepithelial electron-dense deposits 
corresponding to granular IgG staining (1). MN occurs in 
all age groups and races and both sexes worldwide. 
Recently, the association of epigenetic modifications with 
MN has attracted little Interest, but study of these 
relationships may uncover new mechanisms to explain 
the MN epiphenotype. 

Epigenetlcs is the study of inherited modifications in 
gene expression and phenotype caused by mechanisms 
other than changes in the underlying DNA sequence. 



Unlike genetic polymorphisms, non-genetic partial and 
reversible DNA modifications can occur (2). Histone lysine 
methylation, which Is one of the important epigenetic 
modifications, Is believed to be part of a histone code and 
has been implicated In multiple biological processes 
Including gene activation, silencing, X-chromosome Inac- 
tivation, DNA repair, cell cycle control, and DNA methyla- 
tion (3). Methylation occurs at five major lysine residues 
located within the amino-termlnal histone tails (H3K4, 
H3K9, H3K27, H3K36, H4K20), and one lysine residue 
found within the globular histone domain (H3K79). 
Moreover, these lysines can be mono-, dl- or trlmethy- 
lated (4). Among the various histone lysine methylation 
patterns, this study focused on methylation at lysine 9 of 
histone H3 (H3K9), owing to Its association with con- 
densed and repressed chromatin (5). Trimethylated H3K9 
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(H3K9me3) is a posttranslational modification, and is 
commonly linl<ed with both facultative and constitutive 
heterochromatin formation as well as transcriptional 
repression (6). Aberrant alterations in histone lysine 
methylation patterns that change chromatin structure 
could lead to dysregulated gene transcription and disease 
progression. Therefore, it is significant to investigate 
whether there is aberrant H3K9me3 in the nonphysiolo- 
gical MN environment. This may provide important clues 
to assist in the development of new treatments for MN as 
well as to give a deeper understanding of the pathogen- 
esis of MN. 

High-throughput techniques are very efficient meth- 
ods to analyze large-scale gene expression and to 
simultaneously investigate the complex molecular basis 
of pathological processes at the genomic level (7). As 
one of the approaches, ChlP-seq enables the genome 
wide identification of binding sites of transcription 
factors and other DNA-binding proteins (8). With the 
emergence of high-throughput sequencing platforms like 
lllumina Genome Analyzer and SOLiD, and the avail- 
ability of ChlP-grade antibodies, ChlP-seq has become 
one of the most widely used methods for determining 
functional elements in the genome (9). In general, when 
compared with ChlP-chip, ChlP-seq has some advan- 
tages, such as higher signal/noise ratios, lower cost, 
and smaller amounts of ChIP DNA for genome-wide 
analysis (10). 

The understanding of cellular and molecular mechan- 
isms involved in the pathogenesis of MN has come from 
studies carried out in the Heymann nephritis model of MN 
in the rat, which closely resembles the clinical and 
pathologic features of the human disease (11). A number 
of studies have confirmed the relevance of several 
experimental insights to the pathogenesis of human MN, 
but the specific biomarkers of MN have not been fully 
elucidated (7). Hence, in this study, we adopted ChlP-seq 
technology to profile and compare the variations in 
H3K9me3 at the genome-wide level in peripheral blood 
mononuclear cells (PBMCs) from MN patients and healthy 
controls to gain a better understanding of the pathogenic 
mechanisms in MN. 

Material and Methods 

Human subjects 

Twenty subjects were enrolled in the study, which 
included 10 MN patients and 10 healthy volunteers 
(Table 1). The diagnosis of MN was confirmed through 
pathological diagnosis and the presence of other 
clinical evidence. All patients in this study presenting 
to the Nephrology Department of the 181st Hospital had 
renal biopsies positive for MN and were diagnosed to 
be free of any other complications, active infections, 
diabetes mellitus as well as other diseases. Age-, race-, 
and gender-matched healthy controls were recruited by 



advertising. The control group was comprised of 
healthy volunteers diagnosed to be without any inflam- 
matory diseases or conditions. Urine protein, uric acid 
and urea were higher in MN patients, but total protein 
was lower. The Ethics Committee of the 181st Hospital 
approved the study; peripheral blood samples were 
obtained with informed consent from all participating 
individuals. 

Isolation of PBMCs 

Blood samples were obtained from MN patients 
(n = 10) and normal healthy donors (n = 10). Blood 
(5 mL per subject) was diluted with equal volumes of 
phosphate-buffered saline (PBS). An equal volume of 
diluted blood was overlaid on Ficoll-Paque Plus in a 1:1 
ratio and centrifuged at 2700 g for 25 min at 25°C. The 
PBMC layer was harvested and washed twice with PBS to 
remove plasma and Ficoll. We analyzed the total number 
of PBMCs using flow cytometry (Beckman Coulter, USA), 
the final concentration of PBMCs was adjusted to 1 x lO*" 
cells/mL. These samples were then stored at-80°C until 
assayed. 

Chromatin immunoprecipitation 

The Chip was performed according to described 
protocols (Epigentek, USA), with some modifications. 
The preparation of sheared chromatin (cell fixation, 
enzymatic shearing) was performed following the product 
manufacturer's instructions (Epigentek). Briefly, cells 
were fixed with 1% formaldehyde in flasks at room 
temperature for 10 min and subsequently washed with 
ice-cold PBS before a 0.1 -M glycine solution was added 
to stop the fixation. The cells were scraped from the flask 
and homogenized with an ice-cold Dounce homogenizer 
(Active Motif, USA) to release the nuclei. The nuclei were 
then collected and resuspended in nucleus lysis buffer. 
The lysate was incubated with anti-K9 trimethylated 
histone H3 antibody (Epigentek) at 4°C overnight. The 
immune complex was treated with proteinase K for 2 h at 
65°C and DNA was purified with a QIAquick PCR 
purification kit (Qiagen, USA). 

Next-generation sequencing 

All standard protocols for lllumina sequence prepara- 
tion, sequencing, and quality control were provided by 
lllumina (USA). In short, DNA recovered from a conven- 
tional Chip procedure was quantified using the 
QuantiFluor fluorometer (Promega, USA). The DNA 
integrity was verified using the Agilent Bioanalyzer 2100 
(Agilent, USA). The DNA was then processed, including 
end repair, adaptor ligation, and size selection, using an 
lllumina sample prep kit following the manufacturer's 
instructions (lllumina). Final DNA libraries were validated 
and sequenced at 75 bp per sequence read, using an 
lllumina GAIIx sequencer at a depth of approximately 30 
million sequences per sample. 
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Table 1. Demographic and clinical characteristics of the patients with membranous nephropathy (MN) and controls. 





MN patients 


Healthy controls 


Males/females (n) 


6/4 


6/4 


Age (years) 


39.50 ± 15.49 


37.89 ± 13.60 


Total protein (g/L) 


43.06 ± 8.90 


72.85 ± 8.62 


Creatinine (i^M) 


78.33 ± 5.85 


52.01 ± 3.24 


Urine protein (g/24 h) 


1.67 ± 0.86 


0.00 ± 0.00 


Urea (mM) 


5.05 ± 1.95 


2.92 ± 1.17 


Uric acid (iiM) 


375.80 ± 89.48 


121.40 ± 29.37 



Data are reported as means ±SD. 



ChlP-seq analysis 

Chip DNA-end repairing, adaptor ligation, and ampli- 
fication were performed as described earlier (12). 
Fragments of about 100 bp (without linkers) were 
isolated from agarose gels and used for sequencing with 
a Solexa/lllumina 2G genetic analyzer (CD Genomics, 
USA). Sequencing tags were aligned to the bosTau4 Oct 
2007 release of the reference genome using SOAP 2.21 
(BGI, China) and Bowtie, an ultrafast memory-efficient 
short-read aligner (13). We considered those tags that 
aligned uniquely with less than two mismatches. For 
enriched-region (peak) identification (peak calling), we 
used the Model-based Analysis of ChlP-seq (MACS) 
algorithm (14) with the MACS version 1.4.0 software 
(BGI). A simple yet effective technique for the analysis of 
eukaryotes, MACS was designed to identify transcription 
factor binding sites and histone modification-enriched 
regions in ChlP-seq data sets, with or without control 
samples (15). Sequence reads that map to multiple sites 
in the human genome were removed. The MACS 
program was used to determine enriched H3K9me3 
peaks using healthy volunteers as control. The param- 
eters for MACS were: effective genome size = 2.70e-l-09; 
tag size = 49; model fold = 10.30; P value cutoff = 0.05. 
Peaks were displayed through the UCSC Genome 
Browser (http://genome.ucsc.edu/). The enriched peaks 
were annotated with the gene annotation using AmiGO 
1.8. AmiGO is accessible online at the Gene Ontology 
(GO) website (http://www.geneontology.org/) for users 
to obtain the data provided by the GO Consortium (16). 
MEME 4.7.0 was used for motif discovery in the present 
study. 



Results 
ChlP-seq 

Using specific antibodies, we performed a ChlP-seq 
on two samples: histone H3 lysine 9 trimethylation of 
PBMCs from healthy controls and from MN patients. We 
used Bowtie (USA), an ultrafast memory-efficient short 
read aligner (17), to align sequencing reads to the 
reference genome, and the alignment results are pre- 
sented in Table 2. All ChlP-seq data aligned well with the 
reference genome. This process resulted in a total of 
12,899,881 reads and 12,440,646 (96.44%) aligned 
reads, covering 539,358,043 bases. 

Genome-wide distributions of histone H3 lysine 9 
trimethylation of PBMCs from MN patients 

MACS 1.4.0 was used for enriched region (peak) 
identification, with the arguments set as follows: band- 
width =273, model fold = 10.30, P value cutoff =0.05. The 
range for calculating the regional lambda was 5000 bp. 
The ChlP-seq for the MN sample generated 217 enriched 
regions (peaks), total length was 192,047 bp. The peak 
distribution of the ChlP-seq reads of H3K9me3 in MN 
patients exhibited a typical priority distribution pattern 
(Figures 1 and 2), with peak distances ranging from 150 to 
200 bp. 

GO analysis of peaks relative to annotated genes 

To further understand the functions of the annotated 
genes related to peaks, they were functionally classified 
using GO terminology. In terms of the GO database, the 
differentially expressed proteins encoded by these genes 



Table 2. ChlP-seq and alignment results. 



Sample 


Total reads 


Aligned reads 


Aligned reads only 


% Aligned 


% Only 


H3K9me3 


12,899,881 


12,440,646 


11,007,307 


96.44% 


85.33% 



Aligned reads: the number of reads refers to the reference genome; Aligned reads only: the number of reads only refers to the 
reference genome; % Aligned: aligned reads/total reads; % Only: aligned reads only/total reads. 
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Figure 1. ChlP-seq peak distribution (distance) 
of H3K9me3 from membranous nephropathy 
patients. 
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were divided into three groups: cellular components, 
molecular function, and biological processes. The 86 
genes included in the GO analysis are shown in Figure 3. 
From the cellular component categories, we could 
determine the subcellular location of 28 annotated gene 
products that were cellular components; 14 located on 
organelles. Details of the cellular component categories 
are available in Supplementary Table S1. Molecular 
function ontology revealed a binding function in 25 and 
catalytic activity in 12. Genes with enzyme-regulatory 
activity comprised 13.3% of the total. Details of the 
molecular function ontology are available in Supple- 
mentary Table S2. A total of 23 genes were involved in 
cellular processes and 15 involved in biological regulation. 
Moreover, there were 13 genes involved in metabolic 
processes, 10 involved in regulation of biological pro- 
cesses, and 5 involved in component organization. Details 
of the biological process ontology are available in 
Supplementary Table S3. 



Motif analysis of peaks relative to annotated genes 

The detection of functional DNA motifs has been 
greatly facilitated by the availability of high-throughput 
functional genomic data that provide direct or indirect 
evidence forgone regulation (18). MEME is an efficient de 
novo motif discovery tool for large-scale genomic 
sequence data, the motif is believable when E-values 
<0.05. Using MEME 4.7.0, five motifs were searched for 
the ChlP-seq dataset (Table 3, Figure 4). A de novo 
search for H3K9 binding motifs indicated that histone 
modification at different locations changed the sequence- 
based binding preference of histone H3. 

Comparison of H3K9me3 status between MN patients 
and healthy subjects 

By applying the above analysis procedure to the 
sequence results, we found that 108 genes displayed 
significant differences in H3K9me3 in MN patients 
compared with healthy subjects. Among these genes, 
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Figure 2. Genome-wide distribution of peaks 
relative to annotated genes. 
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Figure 3. Distribution of peal<s relative genes by Gene Ontology (GO) analysis (AmiGO 1.8). The lateral axis represents the GO 
terminology. The left vertical axis represents the proportion of the related genes. The right vertical axis represents the number of the 
related genes. 



75 displayed Increased H3K9me3 and 33 genes showed 
decreased H3K9me3 expression. Tlie H3K9me3 altera- 
tions of 15 selected genes that showed the largest 
differences are reported in Table 4. Detailed lists of 
genes with H3K9me3 alterations are provided in 
Supplementary Table S4. Among these genes, DGCR6, 
SNX16, CNTN4, BIRC2 and BIRC3were selected. These 
five genes showed the largest differences, and they may 
be associated with MN. 

Discussion 

Modifications of histone tails were thought to specify 
a code that regulates the expression of genes. Trimeth- 
ylated H3K9 (H3K9me3) is a highly conserved histone 
post-translational modification and is commonly linked to 
both facultative and constitutive heterochromatin forma- 
tion and transcriptional repression. H3K9 methylation was 
initially studied as a marker of epigenetic silencing and 



constitutive heterochromatin, which was found close to 
the centromeres. To cover these regions, we used ChlP- 
seq in this study, which enables the genome-wide 
identification of binding sites of transcription factors 
and other DNA-binding proteins. H3K9me3 was first 
described several years ago (19), but still little is known 
about its subtle interrelationship with other epigen- 
etic modifications and its potential functional signifi- 
cance in human disease. Consequently, we selected 
H3K9me3 as the target of investigation using a ChlP-seq 
strategy, and tested the hypothesis that H3K9me3 
was associated with the pathogenesis of membranous 
nephropathy. 

In the present study, we focused on the trlmethyla- 
tion status of H3K9 in MN patients and healthy subjects. 
The identified candidate genes with significant methyla- 
tion differences are available in Supplementary Table 
S4. Overall, these include genes associated with 
immunity, cell signal transduction, protein transcription 



Table 3. Basic information of motif of H3K9me3 from patients with membranous nephropathy. 



Motif ID 


Width (bp) 


Number 


E-value 


Base sequence 


Coverage 


Motifl 


12 


151 


8.6e-182 


GAVTKGAATGSA 


69.91% 


Motif2 


12 


139 


1.1e-150 


AYRGARTKGAAY 


64.35% 


Motif3 


12 


87 


8.3e-130 


AGAATGCTTCTG 


40.28% 


Motif4 


12 


125 


5.8e-149 


ATYCMMTTSCAK 


57.87% 


Motifs 


12 


71 


7.5e-122 


AAAGAGTGTTTC 


32.87% 


S: C + G; 


R: A-l-G; K: TH-G; Y: C-l-T; M 


AH-C; V: 


A + C + G. 
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^i\T9C^gnaCAG 

ModlS 

Figure 4. ChlP-seq motif logo. The lateral axis indicates locus of 
motif. The total height of the vertical axis reflects the conservation 
of the motif. The height of each base represents probability of the 
base. 

and synthesis of channels and transporters, and extra- 
cellular matrix, etc. 

Among the identified candidate genes, we found that 
DGCR6 displayed increased H3K9me3 in MN patients 
compared with healthy controls. The DiGeorge critical 
region 6 (DGCR6) gene exists as two highly homologous 



copies (DGCR6 and DGCR6L) on chromosome 22q11 
and is deleted in patients with velo-cardio-facial syn- 
drome/Di George syndrome (20). The possible biochem- 
ical activity of the gene product(s) might be gene 
regulation. High levels of the DGCR6 protein have been 
found in ovaries, heart, kidney, adrenal gland, and, to a 
lesser extent, in liver and spleen (20). Sorting nexin 16 
{SNX16) encodes a member of the sorting nexin family, 
which contain a phox (PX) domain. The Phox (FX) 
domain, which is a phosphoinositide-binding domain, 
represents a novel Pl-binding module capable of regulat- 
ing membrane targeting of about 47 mammalian proteins, 
30 of which are tentatively referred to as sorting nexins 
(SNXs). Some SNXs have been implicated in regulating 
membrane trafficking in the endocytic pathway (21). 
SNX16 is a novel 343-amino acid protein consisting of a 
central PX domain followed by a potential coiled-coil 
domain and a C-terminal region (22). Like other sorting 
nexins, SNX16 associates with the membrane via the PX 
domain, which interacts with the phospholipid phosphati- 
dylinositol 3-phosphate. The most documented function of 
the PX domain is to target proteins to membranes 
enriched in Ptdlns(3)P or other Pis (21). The PX domain 
may also participate in protein-protein interactions. Many 
PX domains contain a Pro-rich PXXP motif, suggesting 
that they may mediate interaction with SH3 domains. The 
PX domain may also function to regulate the activity of its 
host protein. Contactin 4 {CNTN4) encodes a member of 
the contactin family of immunoglobulins. In this study, the 
expression of CNTN4 was upregulated. This may be 
associated with disposition of immune-complex involved 
in MN. The presence of glomerular subepithelial immune 
deposits is the distinctive pathologic feature of MN, thus 
supporting the concept of an immunologic origin (23). 
Deposits of IgGI, lgG2, lgG3, and lgG4 were detected in 



Table 4. Selected genes with H3K9me3 alterations between patients with membranous nephropathy and healthy controls identified by 
ChlP-seq. 



Gene ID 


Gene name 


Peak region 


Fold enrichment 


Element 


NM_005675 


DGCR6 


chr22: 18879798-1 8880279 


13.83 


up20k 


NM_022133 


SNX16 


ch r8 : 82754479-82754928 


11.84 


up20k 


NM_175612 


CNTN4 


chr3:3065385-3065914 


10.19 


up20k 


NM_001166 


BIRC2 


chr11:102215301-102215923 


10.00 


up20k 


NIVI_00 11 84940 


C9orf25 


chr9:34459861-34460363 


9.78 


up20k 


NIVI_0 14581 


OBP2B 


chr9: 1 36085386-1 36085951 


9.52 


up20k 


Nl\/I_015210 


KIAA0802 


chrl 8:8708473-870891 8 


8.75 


up20k 


NM_007031 


HSF2BP 


chr21:45081 606-45082287 


7.61 


up20k 


NM_000348 


SRD5A2 


chr2:31810357-31810988 


7.00 


up20k 


NR_027084 


FLJ36000 


chr17:21903894-21905083 


5.99 


up20k 


NR_028492 


CFLP1 


chrl 0:89557757-89558269 


5.31 


up20k 


NIV1_001165 


BIRC3 


chrll : 1 0221 5301 -1 0221 5923 


10.00 


down20k 


NR_030334 


MIR603 


chrl 0:24578544-24579021 


9.78 


down20k 


NM_001 164811 


PET117 


chr20:18130624-18131235 


8.89 


down20k 


NIV1_001 080537 


SNTN 


chr3:63670843-63671536 


6.06 


down20k 
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MN patients. Contactins are axon-associated cell adhe- 
sion molecules that function in neuronal network forma- 
tion and plasticity. The encoded protein is a 
glycosylphosphatidylinositol-anchored neuronal mem- 
brane protein that may play a role in the formation of 
axon connections in the developing nervous system. 
Deletion or mutation of this gene may play a role in 
neurodevelopmental disorders, including autism (24). 
Alternative splicing results in multiple transcript variants. 
CNTN4 has been associated with cerebellar degeneration 
in spinocerebellar ataxia type 16 (25). Disruption of a 
single copy of CNTN4 is sufficient to confer key aspects of 
the 3p deletion syndrome phenotype including develop- 
mental delay (26). Thus, CNTN4 is one of a relatively 
small number of autosomal genes identified as contribut- 
ing to cognitive deficits in a contiguous-gene syndrome. 

We also observed that the other two candidate genes 
BIRC2 and BIRC3 showed a significant decrease in histone 
H3K9me3 in MN patients compared with healthy subjects. 
Baculoviral inhibitor of apoptosis (lAP) repeat containing 2 
{BIRC2) and baculoviral lAP repeat containing 3 {BIRC3) 
encode a member of the lAP family of proteins that inhibit 
apoptosis by binding to tumor necrosis factor receptor 
(TNFR)-associated factors TRAF1 and TRAF2, probably by 
interfering with activation of ICE-like proteases. TRAF2 is 
the main adaptor protein responsible for the signals 
triggered by TNFR2; its degradation exerts a regulatory 
role over the biological activity of the receptor (27). The 
encoded protein inhibits apoptosis induced by serum 
deprivation but does not affect apoptosis resulting from 
exposure to menadione, a potent inducer of free radicals. It 
contains 3 baculovirus lAP repeats and a ring finger 
domain. lAP proteins are negative regulators of cell death 
and play pivotal roles in cellular survival by blocking 
apoptosis, modulating signal transduction, and affecting 
cellular proliferation (28). BIRC2 is essential for maintaining 
endothelial cell survival and blood vessel homeostasis 
during vascular development (29). BIRC2 regulates 
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endothelial cell survival in vivo by modulating TRAF2- 
dependent signaling and, possibly, by ubiquitinating target 
proteins in the TNFR-associated complex, such as NEMO 
and TRAF221. BIRC2 appears to function by modulating 
death-receptor signaling in endothelial cells, possibly by 
regulating the formation and function of TNFR complexes. 
This may contribute to kidney disease progression. BIRC3 
plays an important role in promoting estrogen dependent 
breast cancer cell survival and protecting against TNFct- 
induced cell death (30). The nature of the BIRC3 estrogen 
response element is nearly identical to a consensus 
sequence, and conversion to a perfect consensus 
sequence does not significantly increase estrogen receptor 
activity, suggesting that another factor may be preventing 
estrogen receptor from interacting with this binding site. 

Taken together, we systematically evaluated the 
status of H3K9me3 in PBMCs of MN patients and gained 
new insights into the links between key genes and histone 
methylation in the context of MN. However, as our study 
relies on only 10 patients and 10 controls, there are some 
limitations. Our results may indicate that H3K9 trimethyla- 
tion is involved in nonphysiological MN environment and 
these novel candidate genes may become potential 
biomarkers or future therapeutic targets. Further investi- 
gations with a large sample are needed to clarify and 
validate the roles of identified H3K9me3 candidate genes 
in the pathogenesis of MN. 
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